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Scattering of an electron in quasi-one dimensional quantum wires have many unusual features, not 
found in one, two or three dimensions. In this work we analyze the scattering phase shifts due to an 
impurity in a multi-channel quantum wire with special emphasis on negative slopes in the scattering 
phase shift versus incident energy curves and the Wigner delay time. Although at first sight, the 
large number of scattering matrix elements show phase shifts of different character and nature, it is 
possible to see some pattern and understand these features. The behavior of scattering phase shifts 
in one-dimension can be seen as a special case of these features observed in quasi-one-dimensions. 
The negative slopes can occur at any arbitrary energy and Friedel sum rule is completely violated 
in quasi-one-dimension at any arbitrary energy and any arbitrary regime. This is in contrast to one, 
two or three dimensions where such negative slopes and violation of Friedel sum rule happen only 
at low energy where the incident electron feels the potential very strongly (i.e., there is a very well 
defined regime, the WKB regime, where FSR works very well). There are some novel behavior of 
scattering phase shifts at the critical energies where S'-matrix changes dimension. 



PACS: 73.23.-b, 72.10.-d, 72.10.Bg 
I. INTRODUCTION 



Elastic scattering in one, two and three dimensions is 
well understood In a scattering process there are 
some important physical quantities like scattering ampli- 
tude and scattering phase shift. While the scattering in- 
tensity is directly related to the scattering amplitude, the 
scattering phase shifts are also very important physical 
quantities and the Friedel-sum-rule (FSR) relates them 
to the density of states (DOS). 



FIG. 1. A rectangular quantum biliard or quantum 
dot, weakly coupled to leads. The dotted line is along 
the X-axis and the dashed line is along the y-axis. 



At low temperatures, inelastic collisions are greatly 
suppressed. As a result the phase coherence length of 
an electron can become a few microns. Mesoscopic sys- 
tems are defined as systems in which the phase coherence 
length exceeds the sample size. In such a system elastic 
scattering is the dominant feature and such mesoscopic 
samples can be understood as phase coherent elastic scat- 
terers. With the experimental realizations of mesoscopic 
systems and their possibility of being applied in nano- 
technology and quantum computing, understanding scat- 
tering effects in quasi-one dimensions (QID) has also be- 
come important at present. This is also essential because 
the Landauer conductance formula relates the conduc- 
tance to partial scattering intensities and also one can 
now probe scattering phase shift directly in an experi- 
ment Mesoscopic samples are normally made up 
of metals or semiconductors, and the defects in them are 
generally point defects. Hence we will restrict our anal- 
ysis to delta function potential impurities. 

Recently a new kind of scattering phase shift was dis- 
cussed in QID, in connection with the violation of the 
parity effect jH) . To explain this phase shift and the vio- 
lation of parity effect, here we elaborate some portions of 
Ref. 1^ and explain what are symmetry dictated nodes 
(SDN) and non-symmetry dictated nodes (NSDN) that 
can arise in a QID system. For example let us consider a 
rectangular quantum biliard or dot connected to leads by 
quantum mechanical tunneling as shown in Fig.l. The 
system has reflection symmetry across the x-axis as well 
as the y-axis. Also the x-y components separate and 
spanning nodes (nodes that span across the direction of 
propagation and shown by dashed line in Fig.l) as well 
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as non-spanning nodes (shown by dotted line in Fig.l) 
develop in the geometry, dictated by the reflection sym- 
metries. There are various symmetries that give rise to 
nodes in the wave- function and we call them SDN (for 
example the antisymmetric property of the many body 
wave function result in nodes). But if quantum coher- 
ence extends to some distance inside the leads then one 
can model the phase coherent quantum dot as shown in 
Fig.2. 



leads) 



leads 



FIG. 2. A more realistic model of the quantimr dot of 
Fig.l. 
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(b) 



develop nodes that try to develop across y-direction but 
also act across x-direction and change the phase of the 
wave function in x-direction by tt. We call them NSDN 
because they do not originate due to the symmetry of 
the Hamiltonian. There are many configurations of such 
NSDN 1^, the simplest one was discussed for the stub 
geometry (shown in Fig. 3) in ref. ||^. 

Consider for example two finite one dimensional quan- 
tum wires of equal length, AB and CD placed perpendic- 
ular to each other as shown in Fig. 3 by the solid lines. 
When CD is completely detached from AB as shown in 
Fig. 3(a), then the quantum mechanical wave function in 
AB and CD in the ground state is shown by the dotted 
lines. They are basically the ground state wave function 
in an infinite potential well in one dimension (ID). As 
is known to us, the ground state wave functions are by 
symmetry, even parity states without any nodes, except 
at the boundary. But when CD is attached to AB to give 
a T-shaped stub structure as shown in Fig. 3(b), then CD 
forms a node at C, which is also the midpoint of AB. The 
wave function in this case is again shown by dotted lines 
and the wave function between A and B is no longer an 
even parity state but an odd parity state. The node at C 
between A and B does not originate from the symmetry 
of the Hamiltonian and is not a symmetry dictated node 
(SDN). It is rather forced by the boundary condition in 
the y-direction and is a NSDN. An infinitesimal change 
in the length CD makes this node disappear and then we 
have no node between A and B. The node at C induces a 
phase change by tt and when we join A and B together to 
form a ring-stub system, we also get persistent currents 
without parity effect, as parity of the persistent currents 
is sensitive to the number of nodes in the wave function 



FIG. 3. Two one dimensional quantum wires of equal 
lengths, AB and CD, shown by solid lines, placed along x 
and y directions, respectively. The origin (x=0, y=0) is 
at the mid-point of AB. (a) CD is not connected to AB. 
(b) CD is connected to AB. 



Note that in this case also reflection-symmetry holds 
in the x-direction as well as in the y-direction, but x-y 
components do not separate. So by tuning the bound- 
ary condition in y-direction by a gate voltage one can 



Fano resonances [g| are a very general feature of QID 
P-|lT| systems in the presence of defects and the Fano res- 
onances are characterized by a zero-pole pair of the trans- 
mission amplitude in the complex energy plane. When 
semi infinite leads are attached to A and B in Fig. 3(b) 
then the transmission amplitude of the stub structure 
also has zero-pole pair and one gets Fano resonances . 
At the energy corresponding to the pole a charge gets 
trapped by the scatterer and there is also an energy when 
there is a zero transmission across the scatterer. At the 
energy corresponding to the zero, scattering phase shift 
discontinuously changes by tt due to the NSDN . 

It seems at present that this phase due to NSDN is nec- 
essary to understand the experimental results of Refs. 
[^-0. Initial analysis of the experiments in terms of 
Friedel-sum-rule 113] revealed the shortcomings of ap- 
plying Friedel-sum-rule to the quantum dot. The phase 
change due to NSDN can explain the experimental data 
was first proposed in Ref. [[l3[ and later discussed in Refs. 



The Friedel-sum-rule (FSR) can be stated as |19 



e{E2)-e{Ei)^TrN{E2,Ei). 



(1) 
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In ID, 2D and 3D, the equahty is known to be approxi- 
mate, and is almost exact in the WKB regime where gen- 
erally transport occurs. Here N{E2, Ei) is the variation 
in the number of states in the energy interval [Ei, E2] due 
to the scatterer and [EOl 



(2) 



S being an nxn scattering matrix and e*^», i=l,2,....n are 
the n eigenvalues of the unitary matrix S. In differential 
form the FSR can also be stated as 



8E - 2i8E'<'^'^'^^ 



n{p{E)-p„{E)), (3) 



where {p{E) — po{E)) is the variation of the DOS or the 
difference in the DOS due to the presence of the poten- 
tial. p{E) and po{E) can be found by integrating the 
local DOS p(x,y,z,E) and po{x,y, z, E)) which are re- 
lated to the electron probability. However, if FSR is to 
be useful in mesoscopic systems where scattering phase 
shifts can be accurately measured, and where we are al- 
ways interested to study a finite region of space, then we 
must see to what extent it can give the local DOS. In fact 
in ID, 2D or 3D it does so extremely well and so we ask 
the question that how good it is in QID. We shall show 
that in QID, FSR will fail to give the local DOS as well 
as the global DOS at all energy regimes. 

For a symmetric scatterer in a strictly ID system when 
transmission and reflection amplitudes are denoted by t 
and r, respectively. 



S = 



r t 
t r 



In this case, one can show 



86 _ darg{t) 
dE " dE 



(4) 



which means sum of the phases of the eigen values of S 
is equal to the phase of some particular matrix element 
of S. 

But now we know that in systems that are not strictly 
ID, one can have zero-pole pairs and then Eq.(4) is not 
valid because of the tt phase shifts induced by NSDN. 
Note that for the system in Fig. 3(b), although the scat- 
tering matrix is 2 x 2, the point C is connected to 3 
directions and is hence a QID system. A quantum wire 
with a finite width and only one propagating channel pj| ] 
also is a QID system with a 2 x 2 scattering matrix. For 
such systems, that are not strictly one dimensional but 
has a 2 X 2 scattering matrix 



d9_ dargjt) 
dE ^ dE 



de 



but —^^[p{E)-p,{E)), 



(5) 
(6) 



i.e., when we go from ID (with 2x2 S matrix) to QID 
(also with 2x2 S matrix) then Eq. (3) holds but Eq. 
(4) does not hold. This analysis was presented by Lee 
|16| and by Biittiker and Taniguchi [ |l7| . Their analy- 
sis is restricted to the system in Fig. 3(b) and 5-matrices 
that are 2x2. Eq.(3) is not violated in the presence of 
NSDN and tt phase slips because det[S] — r'^ — and 
if arg{t) changes by tt then [arg(t'^)] changes by 27r or 
and hence det[S] is unaffected by the tt phase slips. 
Thus for such a single propagating channel this phase 
shift is well understood by now. It has been emphasized 
that the multi-channel case also needs to be studied , 
specially since scattering phase shifts can now be probed 
experimentally in the single channel case HI as well as 
in the multi-channel case , but no such study has been 
reported so far. 

In particular, the FSR is very important in condensed 
matter Physics, because from the scattering phase shift 
(that can be determined experimentally) one can know 
the local DOS inside a disordered sample without know- 
ing its internal details. Even theoretically, except in very 
simple situations, the wave function inside the scatterer 
has infinite degrees of freedom and exact calculation of 
local DOS from the exact wave function inside the scat- 
terer may be non-trivial. While a good estimate of local 
DOS from the S-matrix, which is on its own a very use- 
ful quantity, can greatly reduce the complexities. In this 
work we will study the n channel scattering problem in 
a QID quantum wire, with special emphasis on FSR and 
Wigner delay time ||2^,^. Refs. ||l6| and [0 parame- 
terize the 5'-matrix in a particular way (there are in fact 
many different ways of parameterizing the S'-matrix ) in 
which the scattering matrix elements become indepen- 
dent of energy. We will show that this energy depen- 
dence, that are not important in ID play a very crucial 
role in QID multi-channel scattering. Hence in section 
II we will generalize the work of Refs. and for 
real energy dependent 2x2 scattering matrices. The n 
channel case will be analyzed in section HI and IV. In 
section V we will show some novel phase shifts at critical 
energies where S matrix changes dimensions. Section VI 
is devoted to conclusions. 



II. SCATTERING IN ONE-DIMENSION AND 
NEGATIVE VALUES OF DO/DE 

In Fig. 4 we consider a potential that is described in 
details in the figure caption. The quantum mechanical 
wave function or the solution to the Schrodinger equation 
in different regions is also shown and explained in the fig- 
ure and its caption. We will always normalize the incom- 
ing wave- function such that its amplitude is 1. Griffiths 
boundary conditions for this system gives the following 
equations |23 23] (we use 2m — 1 and ?i = 1). 



1 



(7) 
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ae^'=' + be-'''' 



(8) 



where p 



and 



(12) 



ik{l-r) -ik{a-b) = -V{l + r) and (9) 
ik{ae'''' - be-'''^) - ikt = ~V{ae'''' + be-''''). (10) 
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FIG. 4. Two identical delta function potentials sepa- 
rated by a length I. Strength of each potential is V . The 
thick vertical lines denote the positions of the potentials 
and the thin horizontal line is the direction of propaga- 
tion. A plane wave of unit amplitude is incident from 
the left and wave function in different regions (marked 
as I, II and III) is written down in the figure, r and t are 
the reflection and transmission amplitudes, respectively, 
of the entire system and k = ^pE is the incident wave 
vector. The origin of coordinates is shown in the figure. 



{ab* 



2ikx 



ha" 



--2ikx\ 



dx. 



(13) 



Here ^ is a term that arises because of quantum me- 
chanical interference and it can be seen that the inte- 
grand in Eq.(13) oscillates with x. For |^| > 1 (this is 
the WKB regime when the electron does not feel the po- 
tential strongly and is almost entirely transmitted) ^ is 
negligibly small. This is shown in Fig. 5, where we plot 
p (the dashed curve) and p (the dotted curve) . The two 
curves are almost the same for | ^ | > 1 , which means 
Y", being the difference between the dashed and dotted 
curves is vanishingly small above this energy. It is known 
that to get the equality between the LHS and RHS of 
Eq. 11, it is necessary to drop the term ^ |^6|. It is also 
known that this deviation arises because we are consid- 
ering the local DOS rather than the global DOS. If we 
consider the global DOS, i.e., 



p{E) = 



hv 



ip* {x)'ip{x)dx 



where 'ijj{x) is the quantum mechanical wavefunction at 
X, then instead of Eq. (11) we get 



We will first analyze this system in detail and generalize 
the results of Refs. ||l^,|l^ further by considering realis- 
tic energy dependent r and t, that will later help us to 
accentuate the new features that can be observed in a 
multi-channel disordered quantum wire. 

First of all let us calculate the local DOS and global 
DOS to see how much it agrees with dO/dE. Although 
the basic facts discussed in this section is known in the 
Greens function formalism, to the best of our knowledge, 
quantitative disagreement (or agreement) has not been 
shown so far. Using quantum mechanical expression for 
the local DOS integrated over the region II in Fig. 4, i.e.. 



PR = 7- 

nv 



Akx 



be 



-ikx 1 2 



dx. 



it is easy to show from Eq. 3 that (the Eq. below is 
consistent with Ref. jl7| and some hints on its derivation 
is given in Ref [E5|) 



dO 
d{kl) 



Pq 



^ P 



(say). 



(11) 



de 

d{kl) 



Lt 



{W\' + \b\' 



The equality is still approximate although almost exact 
in the WKB regime because the RHS is positive definite, 
while it is well known that the LHS can become negative 
at low energy (non-WKB regime) as will be demonstrated 
below. To get the equality it is necessary to neglect some- 
thing else as we shall soon see. This second thing that 
we shall drop is however not due to the fact that we are 
considering the local DOS. It is an inherent approxima- 
tion of the FSR even when we are considering the global 
DOS and hence also when we are considering the local 
DOS. The pq/l term or the interference term inside the 
scatterer does not arise in the case when Z ^ as in the 
case to be considered in section IV. All the deviation to 
be observed there is due to this inherent weakness of the 
FSR. We shall also see below that this inherent weakness 
of the FSR is negligible in ID, 2D or 3D but becomes 
very formidable in quasi ID. 
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FIG. 5. The system under consideration is shown 
in Fig. 4. The soHd curve gives the exact d9/d{kl), the 
dashed curve gives the p and the dotted curve gives p. 
This plot is done for VP — —5, h — 1, 2m — 1. 



happens at high energy in ID, 2D and 3D) then, 
reduces to the expression — fe-Wi and then therefore, 

|1 — r r^l^ 

= p. It is shown in section III of Ref. Q (see Eq. 6 

and 7 therein), that to relate ^ to the global DOS, one 
has to neglect the energy dependence of the self energy, 
that depends on the coupling of the system to the leads, 
i.e., r' and t' . Thus our results are consistent with that. 
The exact is shown in Fig. 5 by the solid curve. Note 

that in the relevant energy regime (If I > 1), the solid 
curve is very close to the dashed and dotted curves, which 
means FSR works very well for the local DOS as well as 
for the global DOS. But for |f | < 1, deviates from 
p. But since transport effects in weak localization or dif- 
fusive or ballistic regime occur at Fermi energies, that is 
normally higher in semiconductors as well as metals in 
comparison to the energy where the two curves deviate 
substantially from each other, Friedel sum rule is often 
useful in condensed matter to obtain a good estimate of 
local DOS as well as global DOS. 
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FIG. 6. The system under consideration is shown in 
Fig. 4. The plot is oid9/d{kl) versus kl for the system for 
different values oiVl^ . The dotted curve is for = — 2, 
the dashed curve is for Vf' = —2.1, the solid curve is for 
VP = —5, the long dashed curve for VP = —8. We use 
;i = 1,2m = 1. 

One can prove that 



for any energy dependent reflection amplitude r of one 
of the two identical scatterers in Fig. 4, where, r = e*'""'. 
Hence as indicated by Eq. 11, it would be interesting if 
we can obtain a good estimate of p or p' from 6. In the 
appendix I it is explained that if ^ = ^ — > (which 
means the scatterers are non-dispersive and that only 



d6 1 dE is also well known as Wigner delay time . 
In the stationary phase approximation, it gives the time 
spent by the scattered particle at the impurity site. In 
the low energy regime, where dispersion becomes signif- 
icant, the stationary phase approximation is not valid 
and dO/dE can become negative and does not give a 
meaningful particle delay time. In this regime dO/dE 
becomes negative as the phase velocity becomes larger 
than the group velocity and even larger than the veloc- 
ity of light, and although such super-luminous particles 
can be detected experimentally they cannot carry any 
signal or information. In Fig. 6 we show the negative be- 
havior of d6/d{kl). We find that as the strength of the 
impurities is varied, d9/d(kl) can become more or less 
negative (see Fig. 6), maximizing at VP = —2.1 for the 
symmetric delta potentials. The energy regime, where 
dO/d{kl) can be negative remains the same for all V and 
always | f | < 1 . We have checked for all these values of 
V that apart from this insignificant energy range, FSR 
works very well. FSR has a close counterpart in quantum 
mechanics called Levinson's theorem. It is known that 
Levinson's theorem also breaks down in the presence of 
zero energy bound states |Q] that can be degenerate with 
scattering states. 



5 



III. WIGNER DELAY TIME IN 
QUASI-ONE-DIMENSIONS 



Y 



FIG. 7. Here we show a quantum wire of width W. 
The dash-dotted curve is a hne through the middle of 
the quantum wire, and it is also taken to be the x-axis. 
The origin of the coordinates is shown in the figure. A 
delta function potential V{x,y) = jS{x)6{y — Xji) is situ- 
ated at a; = and y — Ui and marked as x . We consider 
scattering effects when the incident electron is from the 
left. The sub-bands on the left of the impurity is denoted 
as 1 for the first mode (i.e., its wave function can be ob- 
tained by putting n=l in Eq.(14) with appropriate sign 
for kn) and 2 for the second mode (i.e., its wave function 
can be obtained by putting n—2 in Eq.(14) with appro- 
priate sign for fc„). Similarly the sub-bands on the right 
of the impurity is denoted as 3 for the first mode (i.e., its 
wave function can be obtained by putting n=l in Eq.(14) 
with appropriate sign for fc„) and 4 for the second mode 
(i.e., its wave function can be obtained by putting n~2 in 
Eq.(14) with appropriate sign for fc„). The impurity at 
X mixes these wave functions to give a scattering matrix 
element from mode m to mode n. 



In Fig. 7 we consider a quasi-one-dimensional quantum 
wire with an attractive impurity at (0, j/i), having elec- 
trons confined along the y-direction but free to move 
along the x-direction. While the states far away from the 
impurity are good momentum states, the impurity can 
mix the different modes and in this region of mode mix- 
ing, the wave function is ip{x,y) = S„c„(x)xti(?/), where 
Xniy) are the transverse wave functions in the absence 
of the impurity and c„(a;) are position dependent coeffi- 
cients that has to be determined by mode matching. The 
confining potential in the y-direction or the transverse 
direction is taken to be hard wall. Thus the transverse 



wave-function is of the form Xn{y) = Sin^{y-\- ^-). For 
a given width W of the quantum wire one can choose the 
energy range of the incident electron such that only two 
modes are propagating, although, all the other modes (in- 
finite in number, showing that the internal wave function 



can have infinite degrees of freedom, which makes it dif- 
ficult to calculate the exact local DOS from the internal 
wave function) will be present but as evanescent modes. 
For example, if the energy of the electron be E then for 
propagation in the n-th transverse mode (in short wc will 
refer this as n-th mode) the wave-function is of the form 



iknX 



(14) 



where k 
1,2,3,-- 



VE-E„, E„ being ^ 
Here we have used h = 



and 
2rn = 



oo. Mere we nave used ti = 2m = 1. To 
have the n-th mode to be propagating it is necessary 
that fc^ > or 



n < — VE. 

TT 



(15) 



Thus we can choose the energy range where there will 
be two propagating modes, i.e., n = 1 and n = 2 satisfy 
condition (15). The rest of the modes (n > ^^/E) will 
be evanescent, whose wave functions are of the form 



Sin—iy + -) 



(16) 



where k„ = y^En — E. These evanescent modes just 
renormalize the scattering matrix elements and drop out 
of the problem. The transmission amplitude from m-th 
incident mode to n-th scattered mode is given by M 



where d — 1 



2d\/kjyi kji 
p p 

- nn ^ ^ -L nn 



(17) 
(18) 



Here denotes the sum over all evanescent modes and 
denotes the sum over all propagating modes. Eq.(17) 
holds only for inter-subband transmission amplitudes and 
all reflection amplitudes. When we say reflection ampli- 
tude we mean the following. For an electron incident 
from the left, all outgoing channels to the left are reflec- 
tion channels. According to this convention tn, ^221 ^21 
and ij^2 a-''^ reflection amplitudes. Inter-subband trans- 
mission amplitude are then obviously and i23- The 
intra-subband transmission amplitudes t^^ and t24 (ac- 
cording to numbering of channels explained in Fig. 7) are 
given by 



t 
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1 + tn 



and to 



(19) 



Here Tnm is the strength of coupling between the n-th 
mode and the m-th mode. If we take the impurity to be 
a delta function potential i.e., V{x,y) = j5{x)S{y — yi), 
and the confining potential in the y-direction to be hard 
wall {V = CO for — 3- > y > and everywhere else 
except the impurity site x ) (see Fig. 7) then 

r„„j = "iSin^{y, + Y^Sin^{y^ + y ). 
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Apart from the two propagating modes we consider 
two evanescent modes and truncate the infinite series 
of evanescent modes (note that although the series is 
strongly converging, the reason for truncating is differ- 
ent, stronger and explained in more detail after Fig. 12) 
in Eq. (18) and so Eq. (18) becomes 

, , Tas . /Fii r22\ , , > , s 

l + ^ + ^+ * Uir + ^ ='^2 say . 20 
2k3 2k4 \2ki 2k2 / 

The lowest evanescent mode (putting n=3 in Eq. (16) 
gives its wave function) has even parity in the transverse 
direction. For a negative impurity potential i.e., 7 < 0, 
it also has a bound state at E = E^h, where E^i, is given 
by the solution of 



1 -\- ^"^"^ -^44 Q 

2k3 2k4 



(21) 



Since E31, < E31, can be degenerate with scatter- 

ing states (the n—1 and n—2 modes are the scattering 
states). The higher evanescent mode (putting n=4 in 
Eq. (16) gives its wave function) has odd parity in the 
transverse direction and this too has a bound state at 
E = Eib, where £^46 is given by the solution of 



1 + ^=0. 
2k;4 



(22) 



Once again depending on 7, E^i, can be degenerate with 
the scattering states. The effect of including more evanes- 
cent modes is just to renormalize the strength of the im- 
purity potential and does not give anything new 
The scattering matrix in this case is 



plane of the paper or out of the plane of the paper. Thus 
among the 16 matrix elements in Eq. (23) we are left 
with only 5 that arc distinct. They are ti 
and ^24- 

From Eq. (17) and (19), 



Ij ^12j ^22; ^13 
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" 2d2fcl' 

_ iri2 

2d2y/kik2 ' 
' *r22 



1 I im _i_ La± _i_ J_22 

^ 2K.3 ^ 2k4 ^ '2fe2 



1-13 



d2 



and 



^24 — 



2k;4 



"2fci 



(24) 
(25) 
(26) 

(27) 

(28) 



Knowing these matrix elements, the scattering matrix is 
completely known and 9 can also be calculated. 

We find some further relationships between the scat- 
tering phase shifts as follows. First of all 



o-rgitii) = arg(t22) = tan 



_i Re{d) 
Im{d) ' 



(29) 



Secondly, when ^ < E31, < i.e., the bound state of 
the 3rd subband lies in the energy range where one can 
have two propagating subbands, then the bound state 
i?3b drastically changes the scattering matrix elements in 
that energy range. So in this energy range ^ to ^ we 
find 



arg{ti2) T 77 = ^ + 



(30) 



S 



^11 '•12 

^21 ''22 



'-SI 



'13 '14 
^•23 ''24 



'32 '33 '34 



'41 '42 '43 '44 



where r2c 



r2c t2c 

he he 



'11 ^12 
'21 '22 



(23) 



Here negative sign is to be taken when E^i, lies in this en- 
ergy range. Otherwise the positive sign has to be taken. 
9 is to be calculated from Eq. (2) using Eq. (23) . Thirdly 
we find 



(31) 



Note that in contrast to Eq. (30) here the choice of ± 
sign is arbitrary. However consistent with this choice is 
the following 



and 



t2c — 



r ' ' 
'13 '1 



tic — tic due to time reversal symmetry and f2c — 7'2c for 
a symmetric scatterer as that considered here. Now once 
again due to micro- reversebility ti2 = ^21- Also ti2 — '14 
because in both t-y2 ^-^id ^14 the density of states in the 
input as well as the output channel is the same, and also 
the incident channel momenta and the outgoing channel 
momenta are the same in the transverse as well as in the 
propagating direction. Also t23 = '41 because transmis- 
sion amplitude should be independent of the position of 
the observer i.e., whether the observer is looking into the 



argitii) ±- = 9- 



(32) 



where once again -I- sign is to be taken when £^3^ is 
present in this energy range and - sign is to be taken 
when absent. 

We thus find very simple analytical expressions for 9 
in the sense that one need not calculate it from a 4 x 4 
scattering matrix but can calculate it from the argument 
of a single matrix element like t^ or ti2 or <22- These re- 
lations are analogous to Eq. (4) in section I obtained for 
purely one dimensional case, i.e., one need not calculate 
9 from 2x2 matrix but one can find it from the argument 
of a single matrix element. 
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FIG. 8. The system under consideration is shown in 
Fig. 7. The plot is of the argument of various trans- 
mission amplitudes {arg{tj^^) in radians) from incident 
channel m to propagating channel n versus EW^^ . The 
solid curve gives arg^tn), the long dashed curve gives 
arg{ti^) and the dotted curve gives arg{t24)- We use 
7 = -10, y, = .21W and = 



In Fig. 8, we plot only the distinct arguments of the 
scattering amplitudes versus energy of the incident elec- 
tron. We find that all of them show negative slopes over 
a very large range of energy and as already discussed, 
such negative slopes give rise to fundamental questions 
in quantum mechanics pT| , p2| . Now in QID we find that 
this negative slope is not restricted to low energy but 
can occur at any arbitrary energy. Notice for example, 
arg(t'i^) and arg(t24^) show larger negative slopes at the 
highest possible energies for two channel propagation. 
The rest of this section will be devoted to understanding 
these negative slopes that at first sight looks very differ- 
ent in nature and character in the three curves in Fig. 8, 
and also to understanding what will happen when there 
are more than two propagating modes. We will address 
the FSR in the next section. 

It is to be noted that among all these scattering ma- 
trix elements tn and exist in the single channel regime 
(i.e., TT^ < EW^ < 47r^) where t^^ is the reflection ampli- 
tude and t]^3 is the transmission amplitude. The phase 
of t'l^ in the single channel regime is known to change 
discontinuously by tt when t'l^ is 0, i.e., has a zero in 
real energy. In the two channel regime if we write from 
simplifying Eq. (27) 



^13 ^ 



fc2(2K3 +g3) +^332 

fe(2K3 +53) + iK3{g2 + agi) ' 



(33) 



where a 



and gs 



Tss\ with 5=1,2,3. then 



ki i"* r44-|-2K4 

interestingly, we see that it has a zero in complex energy 
and not in real energy. 

If we modify the Breit-Wigner line shape formula of 
ID to include complex zeroes and write 



E — Eq 
E — En 



iVr, 



(34) 



where A is a normalization factor, then just as Tp gives 
the scale over which arg[tmbwiE)] increase at E — Ep, Tq 
gives a scale ov6r whicli dTg^tmbw 

{E)] decrease a.t E ~ Eq 
where \tmbw{E)\'^ also shows a minimum at E ^ Eq (but 
not zero). One can check this very easily (let us say, 
when Eo=2, Ep=l and Fq = Fp = 0.5) and so we do not 
demonstrate it here. Now from Eq. (33) we see that at 
an energy which satisfies the condition 



2^3 + 53 = 0, 



(35) 



the real part of the numerator in Eq. (33) is zero. Con- 
dition (35) is the same as the condition (21) for a bound 
state E311 coming from the 3rd subband that is degener- 
ate with scattering states. So, around this energy where 
Eq.(35) is satisfied (lets say at E = E^b = £^0) o-fgii'is) 
will undergo a drop over an energy scale determined by 
the imaginary part, i-e., Fq = K3(72- 

It can be seen in Fig. 9 that |ti3p (dotted curve) shows 
a narrow minimum around an energy EVF^ ~ 84 (which 
is the solution of Eq. (35) or Eq. (21)) and at this 
energy arg{ti^) shows a very sharp drop over a narrow 
energy range determined by K3g2- Hence by decreas- 
ing/increasing this quantity ^3(72 we can make the phase 
drop sharper/broader. 32 can be made smaller in two 
ways, first by decreasing 7 and second by taking the im- 
purity closer to a node in the transverse wave function. 
The plot for a decreased value of 7 is shown in Fig. 10 
and it confirms this. 




FIG. 9. The system under consideration is shown in 
Fig. 7. The solid curve gives arg{t^^) in radians shifted by 
■J radians in the y-direction and the dotted curve gives 

|i]^3p. Both the functions are plotted versus EW'^ using 
Xi — 0, Hi — AhW and 7 = —15. 
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Note that the quantity K3.92 is actually energy depen- 
dent. But in Fig. 9 and Fig. 10 K^g2 is so small that the 
drop occurs over a scale in which K352 is roughly con- 
stant. For larger values of K3g2, the phase drop will be 
determined by a complex competition between K3 and 52 • 
This is shown in Fig.ll. First of all the scale of the phase 
drop becomes so large that any sensitivity to the position 
of the bound state can not be seen. Secondly, Ksg2 can 
not be taken to be a constant over this large scale and 
the enhancement of the negative slope for EW^ > 79 is 
a signature of the fact that here /t3 — > and so K3g2 — » 
as EW"^ increases. 



and demonstarted in Fig. 6, that the energy where the 
negative slopes occur is always for E/V < 1. 



49 54 69 64 



FIG. 10. The system under consideration is shown in 
Fig. 7. The solid curve gives arg(ti^) in radians shifted by 
J radians in the y-direction and the dotted curve gives 

1^13 p. Both the functions are plotted versus EW^ using 



Xi = 0, yi = AbW and 7 = -10 



Similarly if we rewrite Eq.(28) as 

^' ^ ki{2Ks + ga) + iKagi 
~ fci(2K3 + ga) + iK3{gi + 0g2) ' 

where /3 = then it is clear that the behavior of 
arg{t24) will be qualitatively the same. It is indeed found 
in Fig. 8 that the behavior of arg{t2i) is similar to that 
of argit^s). 

6 = ^ln[det[S]] is shown in Fig. 12 as a function of en- 
ergy, for different values of 7. The minimum in 9 follows 
the Esb and so the energy range where the slope of 9 
versus E is negative is determined by the E^i,. Note that 
when E^b goes out of this energy range the 9 versus E 
has a positive slope everywhere. So in Fig. 12, the nega- 
tive slope arises whenever a bound state E^i, is degener- 
ate with the scattering states (n=l and n=2), and non- 
monotonously scatter and disperse the scattering states. 
For weaker impurities in QID, the negative slope occur 
at higher energies and also are steeper as demonstrated 
in Fig. 12. This is in contrast to what happens in ID 



a? 




EW 



FIG. 11. The system under consideration is shown in 

Fig. 7. The solid curve gives arg{ti^) in radians shifted 
by 27r radians in the negative y-direction versus EW^ for 
7 = -47.1371. The dashed curve is for 7 = -25.197. We 
use Xi=Q and j/j = .2\W. 



^ 2.1 




64 69 74 79 



FIG. 12. The system under consideration is shown 
in Fig. 7. The plot is of 6 in radians versus EW^ for dif- 
ferent 7. The dot-dashed curve is shifted by tt radians 
in y-direction for 7 = —47.1371, corresponding E^i, is at 
EW^ = 35 which is less than the propagating thresh- 
old EW'^ ~ 39 of the second transverse mode. The 
dashed curve is for 7 = —25.197, corresponding £^36 is 
at EW"^ — 80. The dotted curve is for 7 = — 15, corre- 
sponding Esb is at EW'^ = 86.606. The solid curve is for 
7 = -10, corresponding is at EVV^ = 87.982. We 
use Hi = .21W and 2;^ = 0. The arrows accentuate the 
positions of the minima that is shifting towards higher 
energies for weaker impurities. 
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We have used two evanescent modes in our calcula- 
tions because one can include as many evanescent modes 
without changing the nature of the negative slopes as 
long as the positions of the bound states E^b and 
remain the same. One can check this that with four 
evanescent modes and 7 = —6.46584, the negative slopes 
are the same as in Fig. 8, which means that the third 
and the fourth evanescent modes just renormalizes 7 
from -6.46584 to -10. The exact renormalization takes 
place according to a formula 7'' = 7"/d, where 7^* is 
the 7 value used here and 7" is the renormalized value 
of 7 when we use m evanescent modes instead of two. 
d = {1 + r^J{2n5) + T2J{2Ke)....r^J{2n,„), where 
r;^„ 7"S'm2^(y, + m-^2_ solving this one can find 
the renormalized value of 7 i.e., 7" that keep the mini- 
mum of any of the curves for the scattering phase shifts 
considered here unchanged. It is worthwhile mentioning 
that at the band edges (i.e., E 39 and 89, in the figures 
considered in this section), the value of any curve is in- 
dependent of the number of evanescent modes, as all the 
modes get decoupled there. In other words, number of 
evanescent modes considered does not change the nature 
of the negative slopes. Only the positions of the bound 
states are important. 

In order to generalize to arbitrary number of propa- 
gating channels we change our notations slightly. For 
electrons incident from the left /right, we call the scat- 
tered channels towards the left/right as reflection chan- 
nels (the rest being transmission channels) and change 
our notation to 

fii = , (36) 

f22 = ^22 1 (37) 

ri2 = i'i2- (38) 

Thus all possible reflection channels are distinct. All the 
intra-subband transmission channels are also distinct and 
they are denoted as in and t22 where 

ill = i'i3 (39) 
and i22 = hi- (40) 

All other scattering matrix elements are equal to one of 
these 5 elements. In this notation the only difference is 
that the lowest channel (n=l) on the transmission side 
is marked 1 instead of 3. So when we say tn we mean 
transmission amplitude from the n=l channel on the left 
to the n=l channel on the right. We find from Eqs. (29), 
(30), (31) and (32) 

^ arg{rrr,n) = ^ (^/n[det[5]]) (41) 

We find the above relation to be true for any number of 
propagating modes. So m and n can take any integer 
value less than or equal to p, where p is the total number 
of propagating modes. For two propagating modes p=2, 
for three propagating modes p=3 and so on. So Eq.(41) 



is analogous to the ID case given in Eq.(4). That is when 
the dimension of the matrix S becomes very large, then 
it is sufficient to consider the argument of a single matrix 
element in order to calculate the complicated quantity on 
the RHS of Eq. 41. In the energy regime where there are 
two propagating channels, the negative slopes in 9 versus 
incident energy curves are determined by E^h, and when 
there are 3 propagating channels then the negative slopes 
are determined by E^h and so on. 

The scattering phase shifts of transmission channels 
i.e. arg{tm,n), where again m and n can take all possi- 
ble integer values less than or equal to p, show sharp or 
gentle phase drops when the scattering states are degen- 
erate with a bound state, depending on the value of the 
imaginary part in the numerator of tmn- lu the single 
channel regime the imaginary part in the numerator is 
zero and phase drops take the limiting value when the 
phase drops are absolutely discontinuous by tt. Just as 
the discontinuous phase drop in single channel case do 
not affect 6 in any way, the phase drops of the arg(tmn) 
also do not affect in any way and behaves similarly 
as arg(frnn)- 



IV. DENSITY OF STATES AND FRIEDEL SUM 
RULE IN QUASI-ONE-DIMENSION 

The local DOS is given by the following expression |^ 

w 

Pr^ dx dy S{E - Ejn.k^) \ ipm,k^{x,y) P 

(42) 

Here E is the incident energy and R is the integra- 
tion region where modes are mixed, m and km are 
the two quantum numbers that define an incident elec- 
tron wavefunction, AmC^'^'^^ Sin^^{y + ^) whose en- 
ergy is Em^km I where we have taken that the electron 
is incident from the left i.e., a: < 0. ipm,km{x,y) is 
the wavefunction in the region of mode mixing and 

^m.k^ix^y) = E„ci"^(^:^")5*«^(y + f )• Here 
Cn^\x,kn) = CnB^'^"^ for Ti — 1 and n — 2 and 
c^\x,kn) = Cne~'^"^ for n > 2; x being greater than 
or equal to 0. The coefficients C„ can be determined by 
using the mode matching technique. The mode matching 
has been done in details by Bagwell Here the delta 
function potential is taken to be extending from — e to 
-l-e which has to be set to be tending to in the end. 
Pqr can be determined by replacing ^pn,k„ {x, y) by the 
plane wave states in absence of the scatterer and doing 
the integration again. 

Thus we find that for any non-zero incident energy 



{p- Po)r = 




10 



T—[\t23\^ + \t24^ + \t25\' 

nv2 



(43) 



Here vi ^ ^,V2 ^ and i,„„ = t„i„ can be 

obtained by solving the matrix Eqs. given in Ref. |^ (see 
Eq. 23 therein and we have used the same notation i.e., 
tmn here is the same as tmn in Eq. (23) of Ref. ||§]). As 
can be seen in Fig. 12 that ^ is negative over a very 
large energy range while (p — jOo)i? as given by Eq.43 is 
positive. 
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FIG. 13. The system under consideration is shown 
in Fig. 7 with 7 = 1. The Fig. shows some important 
scattering probabilities. The solid curve gives and 
it shows that for EW^ > 50, a particle incident in the 
first propagating mode does not feel the scatterer at all, 
and is almost entirely transmitted intra-channel, 
being close to unity. The dotted curve gives 1*22 P and 
once again for EW"^ > 50, it is close to unity signify- 
ing that a particle incident in the second propagating 
channel is almost entirely transmitted intra-channel. So 
EW^ > 50 is the WKB regime where the potential scat- 
ters the incident electron very weakly. The dashed curve 
gives 30 times |?i2 P and shows strong energy dependence 
not only for EW^ < 50 but also around the highest en- 
ergy [EW'^ ~ 89) or in the extreme WKB limit, its ab- 
solute value being extremely small there signifying ex- 
tremely low inter-channel transmission i.e., the incoming 
particle does not feel the scatterer. 



FIG. 14. The system under consideration is shown in 
Fig. 7 with 7 = 1. The solid curve gives Tr{p{E) — po{E)) 
and the dashed curve gives -0.5i-^ln{det[S]). The two 
curves deviate from each other, where ever the curves in 
Fig. 13 are strongly energy dependent. Otherwise they 
agree. 

Note that if we calculate the global DOS by taking the 
integration region to be from —00 to 00 instead of just 
the region R where the modes are mixed then Eq. 43 
remain the same. One will get some extra integrals that 
are indefinite integrals but using the current conservation 
condition it can be analytically proved (see appendix II) 
that they cancel each other. This means the contribu- 
tion to p{E) and that to po{E) coming from outside the 
region R cancel each other in the absence of a term like 
Pq/l as in section II. Thus in this case (the proof is given 
in appendix II) 

piE) ~ poiE) - (p - po)r 

and both of them deviate identically from due to 

strong dispersion, at any arbitrary energy. 

Since the negative slopes are due to the bound states 
supported by the negative delta function potential, one 
may ask what happens for a positive delta function po- 
tential that does not support any bound state. This situ- 
ation is discussed below and it also elaborates the unique- 
ness of the QID, with respect to the violation of Friedel 
sum rule and shows that large violation can occur in the 
extreme WKB limit. Fig 13 shows the energy depen- 
dence of some important scattering probabilities. As in 
ID the scattering probabilities are strongly energy depen- 
dent for EW^ < 50. For EW"^ > 50 all the curves vary 
slowly with energy (non-dispersive scattering). However, 
one scattering probability |ii2p is also very strongly en- 
ergy dependent at the highest energy of EW"^ ~ 89. 
{k^/V » 1). In this regime scattering is almost entirely 
intra-channel as can be seen that l^np and 1*22 P are 
both almost unity. All other scattering probabilities like 
reflection probabilities and inter-channel scattering prob- 
abilities like 1*12 P are extremely small. There are several 
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other scattering matrix elements that are identical to ti2. 
Although being small, they can have strong energy de- 
pendence, or their energy derivative can be very large. 
As can be seen in Fig. 13 that the dashed curve bends 
down with a steep slope at EW^ ~ 89. It is found in Fig. 
14 that there is a strong violation of Fricdcl sum rule in 
the non-WKB regime as in ID and also in the extreme 
WKB regime (EW^ ~ 89 where k^/V » 1 and scatter- 
ing is almost completely intra-channel) quite unlike that 
in ID. It can be seen in Fig. 13 that in the mid-energy 
range, the scattering probabilities are not very energy 
dependent and the Friedel sum rule holds in this mid- 
energy range as can be seen in Fig. 14. But at regimes 
where the scattering leads to strong dispersion, FSR is 
violated. As the strength of the positive delta function 
potential is increased, this regime where the scattering 
probabilities are not strongly energy dependent becomes 
narrower and Friedel sum rule is violated at all energies. 




FIG. 15. The system under consideration is shown in 
the Fig. 7. The plot is of arg{tii) in radians versus EM^^. 
This plot is for 7 = -25.197, x^ ^ and = A5W. 



1.015 



V. PHASE BEHAVIOR AT CRITICAL 
ENERGIES 



'T>> 0.995 



I? 

S 0.985 



Very interesting phase behaviors can be seen at ener- 
gies where the ^-matrix changes dimension. For example 
for E < ^2 there is only one propagating mode and the 

S'- matrix is 2x2. But for E > ^j-, there are two prop- 
agating modes and the S"- matrix is 4x4. The matrix 
element exists on either side of the energy and 
in Fig. 15 we show the behavior of arg{tii) in the energy 
range that includes EW^ = Air^. Note that it exhibits 
a discontinuous phase drop by ^ at EW^ = 47r^. So far 
only discontinuous phase drops of tt has been observed 
but never -j. From the properties of a 2x2 S'-matrix it 
follows that if there is a discontinuous phase change then 
it can only be of tt [IgUij] . So had the S'-matrix been 2x2 
on either side of EW^= 47r^ the phase drop would have 
been TT. But since the S'-matrix is 2x2 only on one side, 
including EW^ = 47r^, i.e., E < the phase drop is 

also one half of tt. li^iP also has a zero at EW^ — 47r^ 
for all possible choice of parameters ||^, and this zero 
is associated with a ^ phase jump instead of a tt phase 
jump. 



30 35 40 45 50 

FIG. 16. The system under consideration is shown 
in the Fig. 7. The solid curve gives |ii3p. The dashed 
curve is after subtracting 5.27183 radians from arg{ti^) 
in radians. We use 7 = +25.197, Xi = Q and yi = A5W. 

Next we take a repulsive S function potential. It is 
known |^ that at critical energies hke EW'^ = 47r^, |i'i3p 
shows discontinuities. Here |ii3p does not have a zero but 
exhibits a discontinuous jump. At these points arg{ti^) 
also shows non-analytic behavior as demonstrated in 
Fig. 16. In this case -^arg{ti^) is discontinuous. 

VI. CONCLUSIONS 

In a multichannel quantum wire with attractive impu- 
rities, negative slopes in the scattering phase shift ver- 
sus incident energy curves can occur at all possible en- 
ergies. For weaker defects it happens at higher energies 
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and the negative slopes are more pronounced. Such neg- 
ative slopes mean super luminescence |^ , ^ that can 
be observed experimentally. Although such a super lu- 
minescent particle will not give any information about 
the particle delay or information delay, they are of inter- 
est because they demonstrate fundamental principles in 
quantum mechanics. Hence Eq. (41) derived in this pa- 
per may be of use to experimentalists and theoreticians. 
Calculations of variation of DOS due to the presence of 
the scatterer show that FSR can be violated at any ar- 
bitrary energy. That is violation of FSR is not restricted 
to non-WKB regime (the regime where anyway transport 
does not occur) as in ID, 2D or 3D. Rather the violation 
is related to strong energy dependence of the scattering 
matrix elements (dispersive behavior). While in ID, 2D 
and 3D the strong energy dependence of scattering oc- 
cur only in the non-WKB regime, in QID there is no 
systematics. A QID system in the extreme WKB limit 
can also exhibit strongly energy dependent scattering and 
there is no definable regime where the FSR will work and 
hence it may not work even in regimes where transport 
occurs. For attractive impurities, weaker the impurity, 
stronger the violation of FSR while for repulsive impuri- 
ties, stronger the impurity, stronger the violation of FSR. 
We also show that the discontinuous phase drops in the 
single channel case have a counterpart in the multichan- 
nel case wherein the drops can be continuous and we pro- 
pose a line shape formula for them in Eq. 34. When there 
is a third channel of escape for the electron, apart from 
the channel along which it is incident and the channel 
where its scattering phase shift is measured, the phase 
drop becomes continuous. However, these phase drops 
do not affect ^lnDet[S] and hence Friedel sum rule. Fi- 
nally, we discuss some novel scattering phase shifts at 
energies where the S matrix changes dimension. 



Here r = e"^, 



and 5*2 = 



kl and k 



T 

T 



2 m 



E. 52 is the scat- 



tering matrix for the free region II of length I between 
the two scatterers. r' & t' are the reflection & transmis- 
sion amplitudes due to one of the two potentials when 
isolated. 

After cascading these three matrices the resultant scat- 
tering matrix of our system becomes 



S = 



1 _ ^12^2 



And so, 



det[S] = 



I _ j,l2^2 



1 - r'^T^ 



{M + N) 



where M = r'^(l - r'V^)^ 



(i) 



(ii) 



(iii) 



and iV = (1 - r'V2) {2t'\'^T^ - t'^T^). (iv) 
From (i). 



ocp \ ocp 0(j) J 



(v) 
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VIII. APPENDIX I 

Considering the symmetric scattering potential in 
Fig. 4, the scattering matrix S of the structure can be 
found by cascading the scattering matrices of different 
parts, i.e.. 



S 



r t 
t r 



Si®S2® S3, 



where Si ^ S3 ^ ^ ^, ^, 



where 



A = r' 



2^ _ j.l2j2 



and B 



I _ j.12^2 

Using (v), (vi) and (vii) we find 
ddet[S] 



A 2r'B 



B 



+2 



+2 



A{1 + Bt 



2r'TAB - B 



2r'^T^B 

\ — J.l2rj-2 



2Tt' 
I _ j,'2^2 



(vi) 



(vn) 



t'2 + r'2t'2r2 
(1 - r'2r2)2 



B 
dt' 



2r'T^B 



dr 



dr' 



(vui) 
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We will neglect the last two terms in comparison to the 
1st term in Eq.(viii) but we will retain the energy depen- 
dence of r' and t' and the 1st term in Eq.(viii) is actually 
varying with energy very strongly. Hence our results cor- 
respond to real potentials and we do not parameterize 
the S matrix in a special way. We apply this result to 
the case of double delta function potential in Fig. 5 and 
illustrate the significance of the last two terms in com- 
parison with the 1st one. We stress that this calculation 
in this appendix holds even if the (5-function potential 
is replaced by the square-well or any arbitrary potential. 
Thus in the regime where ^ ^ and §^ ^ and using 

9t _ „• 



ddet[S] 



= 2 



A ( 2r'B + 



B 



(1 - r'V' 



IT 



At this point we substitute the values of A and B from 
(vi) and (vii) to get 



ddetjS] 
d(j) 



= 2i 



(1 - r'2T2)3 



N 



(ix) 



I 14 -1 

as T = 1 



As 1 - |r'|2 = |t'|2 

Q = 



r'*V*2+ I r' |4| r 



(1+ I r' 



M 
N 



+ 1 



Now substituting the values of M and N from (iii) and 

(iv) 



Q 



-r'*V"(l~r'V^) + |tT(l + kT) 
|t'P(l+|r'|^) 

r'2(l-r-'2T'')-t'^T^(t'^-2r'^) 
-t"r^(t'^-2r'^) 



Using, r'=|r'|e*^'- and t' =\ t' \ e^"* , 
Q= [|rf |tf |r|^e2*^*e2*(^*-^^) - Ir'^if |r| Ve^'^ 



Now, e = —ln{det[S]) 



de de d(j) 



Prom (x), (h) and (ix) 



(x) 



where, £> = (1+ | / | t' |* r^e*'^'- | r' |^ r^e^'^^ 
+2 I t' pl r' |2 r2e2'(^^+^*^ | r' p e^'^O- (xiii) 



Yi det[S] d4> 



Now it follows from unitarity that e^*^^* ^'■) = e"^ = — 1 
and Irl"' = 1 and so 



l-r'V^ f -Kl 



Multiplying the numerator and denominator by 
i_l^/|4 . we get 

0^ 1- I r' |4 (1 - r'2r2)* 1 



1 - r'2r 
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1- 



r' |4 f + 1 



+2|rf l^'l Ve2'(^*+^^) + 2|rY|t'| Ve2'(^*+^^)]/£) 



1- I r 



' |4 



1 - r'V 



where, Q 



/2t-2 12 



(l-r'V^)* 1 



l-|r'|4 f + 1 



(xi) 



(xii) 



(1 - r'*T 



/*2,*2^ ^ 1+ I |4| ^ 14 



1- I r' 



+ 1 



AT ^ 



= [_e2ie* |rf {|tf +2|rf }-(|tf +|rf {|rf +|tf }) 



(g4ie,|^/|2^2 _ e2i(«*+^-)2|rf t2)]/D . 

Now inside the {} brackets if we use the fact that 
\r' \^ + \t' |2= 1, then 



Q = [-e2'^*|rf (l+|rf )-e«^'|i'|V2+2e2'(^*+^'-)|r'|V2]/£>. 



14 



Multiplying numerator and denominator above by 
g-2i(0,+ej) g^^j putting e2''(''' = e*'' = -1, we get 



Q = 



where, £)' = £>e-2^(^'-+^*). 
Again using |r'p + = 1, 



(xiv) 



(xv) 



Now from (xiv) and (xiii) 



i)' = (l + |r'n[|rfe 



'|2Vl7.'|2p-2»et _ |^/|4^2g-2i(9t-e^) 



_|^'|V2g2i(0«-e.)^2|if |r'|V]. 
As = f, 

D' = (l + |r'|2)[|r'|2e-2*«'+r2(|r'|2 + |if )2] 

where of course | r' ^ + | |2= 1. Substituting D' in 
Eq.(xv), 



Q 



^2_ I ^/ |2 ^-2i0r 
^2_|_ I J.I |2 g-2i0t ■ 



Multiplying numerator and denominator of Q by 
and using e^*^^'"^""^ = e"^ = —1, we get from (xi) 

86 1- I r' 1^ 
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IX. APPENDIX II 



The global density of states is given by 
p{E) = f dx f^^dy 

\i^m,km{x,y)f S{E - Em,kJ 



(i) 



where ^m,k^{x,y) = J2„ oh"'\x)xn{y) and Em,k,n is 
the energy of an electron in the leads. 



E, 



+ 



where m = ±1,±2, as 



2mel^2 2me 
there are two propagating modes in the leads. 

As Xn(2/)'s form an orthonormal set, 

/oc 



9 I 



+ 



Here, vi = and V2 = . 

rrie rUe 

Now, 



(ii) 



dxY\ci'\x) = c['\x) 

-OO ^ ' J —OO 



4\x) 



dx 



dx 



+ 



.(1) 



{x) 



dx + ■ 



= Tl(say), 

where electron is incident in the fundamental mode, 

c[^\x) = e*'=i^ + fne-**^!^ for a; < 

= ill 6**^1^ for a; > 
c(^\x)=f 12 e"*''^^ for a; < 



= ii2 e*'^^^ for a; > 



and for n > 2, 



c^rl\x)=tin e'^"" for x<0 
= hn e"''"^ for x>0 

Tl= dx[l + \fn\^ + 2\fn\ Cos{2hx + m)] 

J —OO 

poo pO poo 

+ / dx\in\'^+ / da;|ri2|^+ / da;|ti2|^ 

Jo J —OO Jo 

+ |ii3p + |ii4|^ + 

Here, r/i is defined as fn = |fii|e~"'^ 

Similarly, for electron incident in the first excited mode. 



T2 



/OO . 
-CSO „ 



(a;) 



-OO 

rO 



= f dx[ 

J —OO 



l + |r22p + 2|f22| Cos(2fc2a; + 772) 



POO pO pOO 

+ / (ix|f22p+ / da;|r2i|^+ / dx\i: 
Jo J — 00 Jo 

+ 1^231^ + 1*241^ + , 

Here, 772 is defined as f22 = 1^22 |e~*''^ 



21 



15 



Therefore, 

p{E) = 2 



l + |rii|^ 



hvi 

\h2\' 



dx 



dx 

o 



l + |r22| 
hv2 



dx 



dx 



hvi 

r dx Cosi2k.x + ,.) 
hvi J_oo 

/ dx Cos{2k2X + 

\hl? + \h2? r^^ , |t22p + |t2l| 



hvi 



dx 



+ 7^(|il3p + |il4p + ---) 
+ ^ (1^231^ + |t24|^ + ---) 

Due to time reversal symmetry, fi2 — r2i & ^12 = ^21- 

Inside the square bracket [ ], 

in the 3rd term we put fi2 — f2i, 

in the 4th term we put r2i = ^12, 

in the 7th term we put ti2 = ^21, 

in the 8th term we put ^21 — ti2- 

Therefore, 



l + |ri2p + 1^=22 



hvi 

2_Lli:„.|2 /-O 

dx 



1*221' 



hv2 



dx 



+ 4^/ dx Cos{2kix + f]i) 
hvi J_oo 

+ / dx Cos(2k2X + T]2) 

hV2 J-00 



+ /^(Itl3p+|il4p+---) 



1 

hv2 



{\t23\^ + \t24\' 



Now adding and subtracting the following terms inside 



2|f22| 

hv2 

1 



/ dx Cos{2k2X -\- r]2) 
'J —00 



(|tl3p + |il4p 



' + \t24\' 



Now |rii|2 + |f2iP + |?ii| 



and |ri2p + |f22 



\' + \ii2\' 



|i2ip-l 



l?2iP = l 



Thus, 
P{E) 



2 

hvi 
2\fi 



dx 



2 

hv2 



dx 



/oo 
dx Cos{2kix + 
-00 

/ Cos{2k2X + ri2) 



hv2 

2 

hvi 
2 

hv2 



(1*13 
(t23p 



' + |tl4|' 



l*24|' 



(iii) 



Now 



hvi 



/oo 2 
dx + - — / dx = po{E) i.e. DOS in the 
-00 hV2 J_oo 

absence of scatterer. Also according to some text books 
p?! , for plane wave states the current conservation con- 
dition 



d_ 

dt 



'ip*ip dr 



J.h ds 







can be satisfied if and only if the wave function vanishes 
at ±00 for all energy. Thus the 3rd and 4th term of 
Eq.(iii) of this section is 0. Besides, neglecting these 
oscillatory terms in the leads is very standard in the 
derivation of Friedel sum rule where it is assumed 
that the carrier concentration is so high in the leads 
that the leads are non-polarizable. However, assuming 
non-polarizable leads alone does not eradicate all energy 
dependence of the self energy as thought in Ref. ||l8| (see 
Eqs. 2 and 6 therein, the parameter can be strongly 
energy dependent, which is dispersive behavior). So we 
get 



1*1 



hvi 
we get, 

piE) 



dx. 



1*21 P f° , 1*22 1 



hvi J _ 



dx, 



2 J- 



dx. 



M 

hV2 



dx, 



1 + |2 + 1^21 P + 1*11^ + 1*21 1' 







|2 I Ii;__|2 I I? J2 I lf_ 12 fO 

dx 



dx 



l + |ri2p + |r22p + |?i2|2 + |t2i| 



2k~iil 

hvi 







hv2 

dx Cos{2kix + rji) 



p[E) = p^{E) 



hvi 



(I* 



131 



*24|' 



Thus p(^)-po(£^) - £r(l*13l' 



2 

hvi 
2 

hv2 

= (P- Po)r 



1*14'' 



p +•••) 



l*14| 



+ T-(|*23|^ + *24|^ +•••) 
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